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Abstract 


Fundamentals of classical thermodynamics and chemical species mole balances have been applied to formulate the mathematical equations to 
compute lithium/air cell reversible reaction voltages, heat transfer rate between the cell and its surroundings, theoretical discharge capacity, and 
thermal efficiency. The developed formulation also provides a set of mathematical equations to compute the total cathode volume for a desired 
maximum discharge capacity of a lithium/air cell or battery. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Metal/oxygen or air batteries or cells are unique because the cell cathode active material, oxygen, is not stored in the battery. 
The cell environmental air can supply oxygen to the porous carbon cathode for the half-cell oxygen reduction process. A general 
diagram of a lithium/air cell is shown in Fig. 1. 

Iron, zinc, aluminum, magnesium, calcium, and lithium have been used to develop metal/oxygen batteries. Most of the developed 
metal/oxygen systems employed alkaline electrolytes (e.g., KOH in H2O)). Only the zinc/oxygen battery has become a commercial 
product. Its main use is powering hearing aids. The lithium/oxygen aqueous electrolyte battery suffers from corrosion of the lithium 
electrode because of its chemical reaction with water [1]. A rechargeable lithium/oxygen battery, using a lithium ion conducting 
organic polymer electrolyte, was reported by Abraham and Jiang [2]. They concluded that the overall cell reaction during discharge 
is, 2Li(s) + O2(g) <> Li202(s). In addition to this reaction, Read [3] suggested, based on oxygen consumption measurements, that 
2Liis) + 1/202(g) <> Li2O(s) is another cell reaction that can occur during discharge. Table 1 shows the characteristics of some 
metal/oxygen battery couples. 

From the calculated data given in Table 1, it is apparent that the lithium/oxygen battery, with O2 supplied from the cell’s 
environmental air, has the highest theoretical specific energy of 11,140 Whkg™! of Li. This energy density is higher than that of a 
methanol/air fuel cell, which is 6098 Whkg~! of methanol. 

The lithium/oxygen organic electrolyte battery differs from the aqueous metal/air batteries in that the oxygen reduction product 
is insoluble in the organic polymer electrolytes [2-4]. The peroxide O27~ and oxide O7~ ions are insoluble in an organic electrolyte. 
Therefore, LizO and Li2O are formed at or near the sites where the peroxide ion 02°- or oxide ion O?~ are generated. A summary 
of the experimental work reported in Ref. [3] is as follows. The organic electrolyte formulation was found to have the largest effect on 
the battery discharge capacity and rate capability. The composition of the discharge product, percent of Li2O2(s) or Li2O,s) in the solid 
product mixture of Li2O2(s) and Li2O,s), was dependent on both the discharge rate and the organic electrolyte formulation. The oxygen 
concentration in an electrolyte was dependent on the electrolyte formulation. During the lithium/air cell discharge, the lower oxygen 
concentration in the organic electrolyte resulted in the preferential formation of LizO,,). Measured oxygen solubility and electrolyte 


* Corresponding author. 
E-mail address: sarwan.sandhu @notes.udayton.edu (S.S. Sandhu). 


0378-7753/$ — see front matter © 2007 Elsevier B.V. All rights reserved. 
doi: 10.1016/j.jpowsour.2007.04.006 


S.S. Sandhu et al. / Journal of Power Sources 170 (2007) 196-209 197 


Table 1 
Characteristics of some metal/oxygen battery couples [2] 
Metal/O2 couple Ideal cell reaction Calculated open circuit voltage at 25°C [V] Theoretical specific energy? [Wh kg7!] 
Oz mass included Oz mass excluded 
Li/O2 ALi(g) + Org) <> 2Li20(s)? 2.91 5200 11,140 
Al/O2 4Al(s) + 302g) +> 2Al203¢s) 2.73 4300 8,130 
Ca/O2 2Cas) + Oxg) <> 2CaQ(s) 3.12 2990 4,180 
Zn/O2 2Znís) + Ox) <> 2ZnOs) 1.65 1090 1,350 


* Only the masses of the active reactant materials are included. The specific energy values are given for the cases of inclusion and exclusion of O2. The cell or 
battery mass starts to increase once discharge begins. 
> The reduction of O2 to O° usually occurs in the presence of a catalyst such as cobalt. In the absence of a catalyst the product is the peroxide ion, O7?~. 


Electrical load 


Li Carbon 
Anode Cathode 


Separator 


Fig. 1. Simple sketch of lithium/air system. 


viscosity data were used to determine the oxygen diffusivity in several organic electrolytes [4]. Read et al. [4] concluded that the 
oxygen transport property, diffusivity, is of critical importance to the performance of a lithium/air or oxygen cell. Dolbey et al. [5] 
reported the results of their effort to increase the practical energy density of the lithium/air battery using carbon-air cathodes catalyzed 
by metals. The manganese (Mn) catalyzed cathodes showed the best performance followed by the cobalt catalyzed carbon cathodes. 
The specific discharge capacity of the Mn catalyzed carbon cathodes was determined to be 3137 mAh g7! of carbon at a discharge 
rate of 0.1 mA Citi coin at a cutoff voltage of 1.5 V. Kuboki et al. [6] investigated lithium/air batteries utilizing a hydrophobic ionic 
liquid consisting of 1-alkyl-3-methyl imidazolium cation and perfluoroalkylsulfonylamide anion. 1-Ethyl-3-methyl imidazolium 
bis(trifluoromethylsulfonyl imide) has a higher lithium ion conductivity and prevents hydrolysis of the lithium anode, and showed 
the best electrolyte performance. The lithium/air cell using this electrolyte performed for 56 days in air and showed a discharge 
capacity of 5360 mAh g7! of cathode carbon. 

For the theoretical formulation work presented in the next section, the electrochemical cell is comprised of a lithium foil pressed 
onto a nickel or copper sheet or grid as the cell anode; LiPF¢ in an organic solvent (e.g., PC or a mixture of PC and DME) as 
the electrolyte held in a polymer such as PAN (polyacrylonitrile); and a mixture of carbon grains, binder (PVDF, polyvinylidine 
fluoride), and polymer electrolyte laminated onto an aluminum or nickel grid as the cathode. The organic polymer electrolyte 
membrane separator electronically insulates the cell cathode from the anode. 


2. Cell mathematical formulation 
2.1. Thermodynamic formulation 
The half-cell reactions are given as follows. At the anode during cell discharge: 


iggy = Litin electrolyte) + CGiniti metal) (1) 


198 S.S. Sandhu et al. / Journal of Power Sources 170 (2007) 196-209 
At the cathode during the cell discharge two reactions are possible: 


2Li(inelectrolyte)' + 2€(incarbon) + Org) = LizO025) (2a) 


2Li(in electrolyte)’ + 2€(incarbon) + 502) = Li206) (2b) 


Under the situation of reactions (2a) and (2b) occurring simultaneously in the cathode reaction layer, addition of reactions (2a) and 
(2b) leads to the following overall reaction stoichiometric equation: 


ALi(inclectrolyte)’ + 4€(incarbon) + 3 Or) = Li206) + LizOxs) (3a) 
If reaction (2a) is the only one occurring then the overall cell reaction is: 

2Liys) + Orgy Z Li20xs) (3b) 
and if reaction (2b) is the only one occurring then the overall cell reaction is: 

2Liis) + 40x) = Li20%) (3c) 
and if both reactions (2a) and (2b) occur then the overall cell reaction is: 

4Liw + 3029) = Li20%) + Li2029 (3d) 


The thermodynamic fundamental criterion of reaction equilibrium is applied to reactions (1) and (2a) to obtain the mathematical 
equations for the prediction of the reversible half-cell reactions that occur at the cell anode and cathode during the cell discharge. 
The condition of equilibrium for reaction (1) is: 


chem __ _,echem echem,A 
LLi(s) = HY i+ (nelectrolyte) + Me-(inLi) (4) 


The condition of equilibrium for reaction (2a) is: 


echem echem,C chem __ ,,chem 
ZUL i+ (inelectrolyte) + He~ (n carbon) + #Ox(g) = HLiz02(8) (Sa) 


chem chem che: 
where Hiie) > KOE) and Hiis 


electrochemical potential of Li* ion in the electrolyte, and u 


m . . . . . echem . 
O(s) are the chemical potentials of Liis), O2(g), and Li2O/s), respectively. HT it Gn electrolyte) is the 


echem, A d echem,C 
e7 (in Li) e7 (in carbon) 
chem 


electron in the anode and cathode reaction layer, respectively. For a charge neutral species (i), p; 


are the electrochemical potentials of an 


is expressed as: 
ye = u? + RT In(a;) (5b) 


where pe is the standard-state chemical potential of a species (i) at the cell temperature and a; is the activity of species (i). The 
electrochemical potential of Li* in the electrolyte solution is expressed as: 


echem — , chem S_ ,,0 S 
T+ Gnelectrolyte) — HLit (inelectrolytey + Lit TO” = MLitGnelectrolyte) + ÈT In(4y j+(inetectrolyte)) + Zit FO (5c) 
where ua ünelectrólyie) is the standard-state chemical potential of Li* ion in the cell electrolyte solution, AL it (in electrolyte) the activity 


of Li* in the electrolyte solution at the cell electrolyte solution pressure and composition conditions, zņ;+ the number of positive 
charges on the lithium ion and has a value of +1 and @° is the electric potential in the electrolyte solution just outside of the outer 
Helmholtz plane (OHP) of the double charge layer between the bulk electrolyte solution and an electrode. The electrochemical 


; hem, A hem,C F 
potentials, u°% (inti) and uec E arbon) ate expressed as: 
hem, A h A 0 A 
Me Gn Li = Me-Gn Li) + Ze- FO” = Me-(inLi) + RT In(4e-(inLi)) + Ze- Fo (6a) 
and 
echem,C _. chem,C C_ 0 Cc 
e~(incarbon) ~ /e~(in carbon) t Ze- Fo“ = Me-(in carbon) + RT In@e~(in carbon)) + ze- Fob (6b) 


where ie GnLi and u? are the standard-state chemical potentials of an electron in the lithium and carbon electrodes, 
i in Li) e~ (in carbon) TO . ane 

respectively, at the cell temperature; ae-(inLi) ANd Ge-(incarbon) are the electron activities in the lithium and carbon electrodes, 

respectively; Ze- is the number of charges on an electron and has a value of —1; and ¢^ and #° are the electric potentials in the bulk 

regions of the lithium anode and carbon cathode, respectively. The electron activities, ag-(in Li) ANd Ge-(in carbon)» are usually set equal 
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to 1. The information provided in Eqs. (5b), (5c), (6a), and (6b) is used in Eqs. (4) and (5a) to obtain, after algebraic simplification, 
the expressions for reversible reaction electrode electric potentials as follows: 


0 0 0 
(gh E 45) = pA = PL it Gn electrolyte) + He- 7T MT i(s) rs (+) In (e eroina ) (7a) 
F F aLi(s) 
papel ( RT ) 7 (“acest | (7b) 
F AL i(s) 
where E^? is the standard-state electric potential of the anode reaction (1). 
p 
0 0 0 0 
(¢° — g5) = ES = (2 + 2H 1+ dn electrolyte) + He- — rons ) Ge In ( ae (8a) 
= = 2 
2F 2F Oxe) AL it (in electrolyte) 
EE = EC? _ (=) In ( AL i707) (8b) 
= 2 
2F 409g) Fy + (in electrolyte) 


where E©° is the standard-state electric potential of the cathode reaction in (2a)with LizOx(s) considered to be the only discharge 
product formed. 

If a pure solid phase film of LizO2 is assumed to be formed on the carbon cathode aLi Ox is set equal to 1. Furthermore, ue, 
the standard-state chemical potential of an electron is also set equal to 0. Also, agi(s) is set equal to 1 since Lis) is assumed to be a 
pure solid phase. The resultant equations to compute reversible anode and cathode reaction electric potentials are: 


RT 
E^ =o" — 95) = E^? + (=) In(4_ j+ (in electrolyte)) (9) 


and 


RT 1 
EC = (4° — #5) = EW (SF) m( 7 ) (10) 
40g)" + (in electrolyte) 


By subtracting Eq. (9) from (10) the reversible electric potential for the overall cell reaction given in Eq. (3b), where Li2O2çs) is the 
discharge product, is obtained. 


RT 
— pC pA (aC AS) aA ASCO pA) —1 
E = E- — E“ = ($~ — °) — ($^ — °) = (E E™") (SF) mas) (11a) 
RT 
= (EC? — FA% + (3) Indore) (11b) 
RT 
= E? + | — ) Inao) (11c) 
2F g) 
where 
1 
0 C,0 A,0 0 0 0 
ee ee (zF) (Hong + 2HLig) T PLi) (hid) 


is the standard-state electric potential of the overall cell reaction (3b). The activity of a gaseous species (i) in a gas phase mixture is 
given by: 


TE iad (11e) 

l l po 
where ĝ; is the fugacity coefficient of a species i in the gaseous mixture to take account of the energetic interactions between the 
molecules resulting in nonideal gas behavior. For gas phase mixtures at pressures less than 10 atm, it is quite reasonable to assume 


that @; = 1. Finally, the reversible electric potential for the overall cell reaction (3b) is: 


E= E+ liin ) (12a) 
= OF Ox2) 
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RT P 
E=E? +n (22: ) (12b) 


2F poe 


which gives the reversible electric potential for the lithium/O2 cell assuming reaction (3b). 

Using the above procedure, one can obtain the following equations to predict the reversible electrical potential of the overall 
cathodic reaction (3a) and that of the overall cell reaction (3d) for the case of formation of Li2O2(s) and LizOQ,s) simultaneously at 
the carbon cathode during cell discharge: 


RT ALi Ors AL in Ores 
40 S S 
EE = E — electrolyte = ESS = ( ) In ( a i20) “Lig 0 ) (13a) 


4F AL i+ Gn electrolyte) O2) 

C0 RT 4 L5 

E simul F ( 4F ) In(@y i+ (in electrolyte) Oze) ) (13b) 
since aLi, Og) = ALin02) = 1 for the assumption of the formation of pure solid phases Liz Os) and Li2O,s) and where the standard-state 
electric potential of the overall cathodic reaction (3a), Ex is given by: 

0 0 _ 0 _ 0. 
C,0 Matra electrolyte) + LSO FL in0) ML in Or) 
BS = (13c) 


simul —— 4F 
Now the reversible electric potential, Esimu, for the overall cell reaction (3d) is given by (assuming Li2Oçs) and Li2O2çs) form solid 
crystal grains): 


RT ALi2 Og) ALi2Or(s 
Esimul = (HF — bs) — (^ — bs) = E? mu — (FF) In (mages | (14a) 
7028) 
RT 
= Fo nu + (FF) In(ag3,)) (with aLiro4 = aLi = 1) (14b) 


where ae is the standard-state electric potential of reaction (3d) and is given by: 


0 0 0 0 
Eo = 4u Lig t 1.5u Org — FLi0w) 7 FLO% 
simul — 4F 
For the situation of a homogeneous mixture of the Li2O2(s) and Li2Q,s) in the solid phase in the cathode reaction layer, the activities 
AL i706) and AL i706) should be computed from: 


(14c) 


AL i706) = YLi20%)XLi20%) (15a) 
ALi20x49) = VLiz02(5)XLi202¢5 (15b) 


where VLi20) and YLigOx) are the energetic interaction coefficients to account for the nonideal behavior of the solid phase mixture 
of Li2O(s) and Li2O2çs) and where x,j,09, is the mole fraction of Li2O2(s) and x; ;,0 is the mole fraction of LizO,.). Eq. (14a) coupled 
with Eqs. (15a) and (15b) provides the cell method to determine the solid phase energetic coefficients experimentally. The activity 


of the Li* ion in an electrolyte solution, ay ;+¢in electrolyte)» 18 expressed as follows: 
ay jt = YLit mit (16a) 


where my ;+ is the molality of the Li* ions (g mol of Li* ion kg~! of solvent), and YLit.m İs the activity coefficient of the Li* ion to 


account for the electrolyte nonideal behavior due to its interaction with the surrounding charged and uncharged particles. If yg i+ m is 


-1 


considered to be a dimensionless quantity; then, my ;+ must be divided by mọ = 1 g-molkg™ solvent. Then, a; ;+ may be expressed 


as: 
My j+ 
ay it = YLit,m 3 (16b) 
mo 
If the Li* ion concentration is given in molar units, it can be expressed in molality units using the following expression: 
Chit 
m+ = H (17) 
PO 


where Cy, ;+is the lithium ion molar concentration, g mol of Li* ion per liter of solution, and pọ is the density of the pure solvent at 
the solution temperature and pressure conditions. The volume of the liquid phase solution is usually assumed to be equal to the pure 
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solvent volume. Analytical expressions to compute activity coefficients of individual ionic species such as Li* ion and mean molal 
activity coefficients of electrolytic solutions can found in Refs. [7,8]. 
The variation of the reversible voltage of the overall cell reaction (3b) or (3d) is given by [7]: 


OE AS 
— |} = (18a) 
where AS and n are the entropy change and g mol of electrons involved per g mol occurrence of the overall cell reaction (3b) or (3d). 


Over a small temperature range, the assumption of constant entropy change of a reaction is usually justified, and integration of Eq. 
(18a) with respect to temperature leads to: 


AS 
E=Eij= C=) (18b) 


where E; is the cell reversible reaction voltage at temperature Tı and £ is the reversible voltage at temperature T. If the entropy 
change varies significantly with temperature, then an integration of the temperature dependent entropy function in Eq. (18a) should 
be carried out. The reversible voltage change for the overall cell reaction (3b) or (3d) with pressure is given by: 


E — E1 = 


2.303 ANRT P 
lo o( 2) (18c) 


nF Pi 


where AN is the change in g mol of the gaseous species [7] due to the 1 g mol occurrence of an overall cell reaction in the forward 
direction where n g mol of electrons are involved in each half-cell reaction at each electrode. Note that AN is —1 and n=2 for 
reaction (3b) and AN is —1.5 and n=4 for reaction (3d). 

The thermodynamic equilibrium constant, K, for an overall cell reaction (3b) or (3d) is given by: 


nFE° 
K = exp RT (19a) 


where E° is the standard-state potential of the cell reaction at the cell temperature, T [K]. K corresponds to the point at which the 
driving force for the further electrochemical reaction is zero. The equilibrium constant, K, is related to the activities of chemical 
species taking part in the overall cell reaction at the cell equilibrium state. That is: 


K= (Te) (19b) 


where v; is the stoichiometric coefficient and a; is the activity of a chemical species taking part in the overall cell reaction. For 
reaction (3b) ULi«) = —2, VOxg) = =l, ULin0x) = L, aig = aLi) = 1, and A072) = YO? P/ P? (for P < 10 atm) and for reaction 


(3d), ULigs) = —4, VOx) = —3/2, VLin 075) = 1, VLi2 0) = 1, aLi) = AL in Oc) = AL in 0x5) = 1, and 10 9/9) = YO> PJ P? (for P< 10 atm). 
Using Eq. (19a), one can calculate the reaction extent to which a reaction can proceed theoretically for a given temperature, pressure, 
and internal reactant composition conditions. 


2.2. The species mole balances and heat transfer 


A lithium/air cell is assumed to be operating under isothermal, i.e., constant temperature, and isobaric conditions to deliver total 
current /* [A] at the cell voltage E°"! [V]. The oxygen supply for the lithium/air cell is assumed to come from the surrounding air at 
temperature, Tsur [K], and pressure, Psy;- [atm]. In the thermodynamic language, the cell is labeled here as the system to be analyzed 
as shown in Fig. 2. It is assumed that there is no mass loss due to solvent evaporation and thermal degradation of components. 


Cell 


(system) O2 from dry air 


Fig. 2. Schematic representation of the overall lithium/air system. 


202 S.S. Sandhu et al. / Journal of Power Sources 170 (2007) 196-209 


2.2.1. The chemical species mole balances 
During the cell discharge period, the rate of increase of the total cell mass is given by: 
diss 
dt 
din\y. 


dt 


= (rate of oxygen mass consumed per unit time in an overall chemical reaction) (20a) 


= (rate of oxygen mass supplied per unit time from the cell surrounding air) (20b) 


Now if reaction (3b) is the overall cell reaction then: 


dmt It 
sS = Ga) (21a) 


and if reaction (3d) is the overall cell reaction then: 


Oi. 3 ri (21b) 
d \8/F 


Masses of the electrolyte, carbon, binding material, and gases in the dissolved state in the electrolyte are assumed to be invariant 
with respect to time. Furthermore, change in the amount of oxygen, and nitrogen present in the electrolyte solution in pore volume 
of the porous carbon cathode is assumed negligible. 

If the overall cell reaction is as given in reaction (3b), then the rate of change of g-mol of Liçs) in the anode and Li202ç(s) in the 
cathode are given by: 


dN. ; It It 
lo 2( ) z (22a) 


dt 2F F 


dM. Tt 
Li20 xs) oe (22b) 
dt 2F 
On integration of Eqs. (22a) and (22b) with respect to time one obtains equations to determine the amount of lithium and lithium 


peroxide in the cell at any time, ¢, for cell reaction (3b): 


It 
Nig = Nigga = F‘ (23a) 
It 
Niir = NLi2029),0 + zF" (23b) 


Similarly, for cell reaction (3d) the amounts of lithium in the cell anode and lithium peroxide and oxide in the cathode, for the 
situation of equimolar formation of Li20O2(s) and Li2Oç¢s), are given by: 


It 
Nig = Migo — p (24a) 
It 
NL = NLi09).0 T JF’ (24b) 
It 
NLirOow = NLO% + IF‘ (24c) 


For the case of different molar formation rates of Li2O2(s) and Li2Oçs), the following approach can be used in the absence of 
experimentally obtained intrinsic kinetic rate laws for their formation rates. Here the total current, J‘, is the sum of the currents, 
Tno% and Trog’ associated with the formation rates, Rino% and Rind? of Li2O2(s) and Li2Q(s), respectively. That is: 


I= Tino% + Tino) (24d) 
Also, 


t t t 
irog = CFRE Ow = Rin dw _ 


I ~ (2F)R! ~ RI Gee) 
Li20xs) QF) Li202(s) Li202(s) 
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where r is the ratio of the solid product species molar formation rates that can be adjusted. Solution of Eqs. (24d) and (24e) for 
Tiro and Kino leads to: 


It 
t = 
ERON = 1+ r (24f) 
t r t 
Linde = (, i -) I (24g) 


The amounts of lithium peroxide and oxide in the cell cathode, for the case of their different molar formation rates, are given at any 
time, t, by: 


1 It 
NL) = NLi202%9),0 + Ir (5) t (24h) 
r It f 
NLirow = NLO% + Tar \(2F t (24i) 


Eq. (24a) is still valid for the amount of solid lithium in the cell anode. Eqs. (24h) and (241) reduce to Eqs. (24b) and (24c) for r= 1; 
the case of equimolar formation rates of the solid product species. 


2.2.2. Cell heat transfer 
The following general enthalpy balance equation for constant temperature and pressure conditions is applied over the lithium/air 
cell as the system: 


[rate of change in the system thermal energy content due to change in the amounts of the solid phase reactant and product species 
at constant temperature and pressure] = [rate of thermal energy production because of occurence of the overall cell reaction] 
+ [rate at which enthalpy enters with the gaseous reactant oxygen] 
+ [rate of heat exchange between the system and its surroundings] 
+ [rate of work exchange between the system and its surroundings in the form of electrical energy] 
+I 


rate of work exchange between the system and its surroundings due to the system volume change at constant pressure] 


(25) 

and this results in the following equation for the reaction (3b): 

d I' I ie has 

a |e MiCea(T — To) | = | 5p C ABr) } + | Hoson, zp | +O + We + Ways. (26) 

i 
Under constant temperature condition, Eq. (26) becomes: 
I I I' ; ; 
(5) (CELO — 2CP.Lig (T — To) = (<7 (-Atnn)) + [Hon az] +Q + We + Wsys. (27a) 


where C P.LinOn(s) and C P.Lig are the solid phase heat capacities of Li2O2çs) and lithium, respectively, on a per g mol basis; T and To 
are the system and reference temperatures, respectively; Q the heat exchange rate; W, the electric work exchange rate; and Was, 
is the rate of exchange of work between the system and surroundings due to the net change in the system volume at the constant 
pressure. Derivation of the left-hand side of Eq. (27a) from the left-hand side of Eq. (26) is as follows. Under constant temperature 
condition: 


dNLi,) AN Li. 024s) C 


dt 


d d 
dt Gy LN Lies CP Lig + Nin 095) CPLinOx) (T = To) = | CP Lig + P,Li202%) 


XON:Cr:T — m) 
It It 
x (T — To) = \(-3) CP Lig + (5) Crund.| (T — To) 


It 
~ (3) (CPLix0x5) — 2CP.Lig (T — To) le 
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For reaction (3b) the enthalpy of reaction is: 
AFrxn = (Htr — 2 ALi, — Horo) (28) 


It should be noted that the quantities Q, We, and Weys. have positive number values if heat and work enter the system and negative 
values if heat and work leave the system. 

During the period of cell discharge delivering total current at the level J‘ [A] at the cell voltage E 
given by: 


cell 


actual [V] the electrical work is 


W, = — IES (29a) 


actual 


It is assumed here that the discharge product, LizO(,), due to the overall cell reaction in Eq. (3b), is formed in the pores of the 
carbon cathode layer. The cell stops delivering electric power once the cathode pores are filled with the discharge product or even 
before the filling of pores with the solid product if oxygen is unable to reach the carbon-electrolyte interface active sites for the net 
cathodic electrochemical reaction to occur due to pore mouth plugging with the discharge product. The system volume change is 
due only to the lithium metal consumption and Li2O2s) formation. Therefore: 


Wsys. = (system volume change per g mol of reaction (2b) occurring)(g mol of reaction (2b) occurring per unit time) Psys. 


M 2Mu\ / I 
= =O = ( ) Bi (29b) 
PLiOx)  PLig 2F 


Substituting for We and Ways, from Eqs. (29a) and (29b), respectively, into Eq. (27a) and further simplifying the algebraic expression 
leads to: 


: It 
Q = (=) AF xn — Hozo T, + 2 FES% + 


Eq. (30) is the expression to calculate the amount of heat required to be removed for the isothermal operation of the cell at a 
temperature, T, during its discharge. 
Enthalpy of a species at any temperature can be computed as follows: 


Mrio  2Mii 


Psys. + (CP, Liz029) — 2CP.Li¢ MT — To) (30) 
PLi2 O26) PLi 


F 
Ais) = Hiis), To + | Cpos),i dT (31a) 
To 
T 
= AH) em + | Cp(s),i dT (31b) 
To 
T 
Hig) = Hig), To + | Crcg),i dT (31c) 
To 
T 
= AHigg)¢.7) + | C Peg),i IT (31d) 
To 
where AH? 


i(s),f,Ty and A Ho), n are the standard-state enthalpies of formation of the solid state and gas phase species from their 
respective elements at the reference temperature of To = 298.15 K and Cp,;’s are the chemical species heat capacities. Following the 
above procedure one can obtain the following expression to calculate heat transfer between the cell system and its surroundings when 
the overall cell reaction is that given in Eq. (3d), where Li2Oç¢s) and Li2zO2,s) are both formed during the period of cell discharge. 


. It MLio Mro 4Mii 
Q= (Fe) [aera Stosor +4 FE Sinai t ( aioe | Psys. HCP Lino +CP.Liz0q —4CP Lio T- To 


4F PLO) PLO% Plig 
(32a) 
where the enthalpy of reaction (3d) is given by: 
AHrxn = (Hiro + HLirow — 4Hig — 1-5 HOx.)) (32b) 


In the derivation of Eq. (32a), it has been assumed that the discharge products Li2O¢(s) and Li2O2çs) are formed in the pores of the 
carbon cathode. When the pore mouths in the carbon cathode are plugged and the reactant oxygen is prevented from reaching the 
active carbon-electrolyte interface sites or when the pores are completely filled with the solid discharge product species, the cell 
stops delivering electrical power. 
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cell 


Actual cell voltage, Æ$ftua When the cell is delivering total current I‘ [A] to an external load, is given by [9]: 


Ea = Erev — (In| F nS] T |A@ionic| + | Agelectronic|) (33) 


where Eey is the reversible cell voltage either for cell reaction (3b) or for (3d). Eq. (33) is strictly valid for the low current density 
condition where resistance to mass transfer of a reactant, for example, oxygen, to the active sites at the cathode-electrolyte interface is 
negligible. At high current densities, resistance to mass transfer of a reactant species to the reaction active sites would also influence 
the cell performance. Then, the cell voltage loss associated with the reactant mass transfer should also be accounted for as explained 
in the reference [9]. Eqs. (12b) and (14b) can be used to compute Erev for reaction (3b) or (3d), respectively. The overvoltages at 
the anode and cathode, |7“| and |7©|, respectively, should either be determined experimentally or from the Butler-Volmer kinetic 
equation [8] written as: 


_ (1 = B)F in| —BF\n*| 
= 4 
i= iĝ (ex ( RT exp RT (34a) 
which is for the lithium oxidation process at the anode, and: 
c C 
PET bnF\n-| —( — n F|n™| 
= 4b 
lil = 151 (cx ( ar) Tex RT (34) 


which is for the oxygen reduction process at the cathode. Here, the symbols i and |i| on the left-hand sides of Eqs. (34a) and (34b) 
represent the net currents per unit electrode-electrolyte interfacial area at the cell anode and cathode electrode, respectively; |n‘4|, 
In@|, |i C] represent their respective absolute or positive number values; 6 = 0.5, n = the number of electrons involved in an elementary 
reaction step that controls the overall oxygen reduction process at the cathode. Assuming the oxygen reduction process follows the 
reduction mechanism scheme in the presence of Li* ions similar to that in the presence of H* ions, one may assume n= 1 in Eq. 
(34b). in and i are the exchange current densities at the lithium anode and carbon cathode, respectively. The voltage loss due to 
the ohmic resistance to Li* ion movement under electric field effect in the electrolyte solution is given by: 


i eom/electrol te 
|Agionic] = === (35) 
Kelectrolyte 
where igeom = I t/ Ageom [A Gti Ageom is the geometric area of the electrolyte separator, electrolyte the electrolyte path length 


[cm], and electrolyte the electrolyte conductivity [S cm7!]. | A@electronic| the voltage drop due to electronic migration in the lithium 
metal and carbon cathode, would be small relative to | Adionic| and |7©| and may be neglected. 


2.3. Theoretical discharge capacity and thermal efficiency 


Theoretical discharge per unit mass of the active material, lithium, is: 


F 
(theoretical = Mu (36) 
where My; is the atomic weight of lithium, 6.939 g g mol~!. Note: Eq. (36) applies to both reactions (3b) and (3d). 
The theoretical or maximum thermal efficiency under cell reversible operating conditions is given by: 
AG 
Nthermal, max = AH (37) 


where AG = )>,ujG; and AH = X ;v;G; are the reaction Gibbs free energy change and the enthalpy change, respectively, per 
g mol of reaction (3b) or (3d) occurring at the cell temperature. 


2.4. Cathode volume estimation 


In this section, the mathematical expressions to compute volumes of the various components and the total volume of the composite 
carbon cathode are provided. It is assumed that the porous cathode is composed of carbon grains; a binder, such as PVDF; and a 
polymer electrolyte, such as LiPF¢ in a mixture of propylene carbonate and 1,2-dimethoxyethane held in PAN (polyacrylonitrile). 
Maximum cell discharge is allowed to the extent that the discharge product, solid lithium peroxide or a mixture of solid lithium 
peroxide and lithium oxide, resides in the cathode pore volume. Also, it is assumed that a small pore volume would be available for 
air or oxygen access to the cathode reaction sites even by the end of a desired maximum discharge from the cell. 
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If the desired discharge maximum is qj“, the amount of lithium consumed is: 


—AN i = (38) 


where q'j""* [C], and is the maximum discharge desired from the battery. The decrease in the lithium metal volume in the anode is 


given by: 
(-AM My qg” Mii 


-AV = = (39) 
Ho pLi FOLi 


where pr; is the density of solid lithium metal [g cm73]. 
If Li2O2z) is the only solid discharge product, as described by Eq. (3b), being formed in the porous composite cathode, then the 
amount of Li2O2(s) deposited is given by: 


—ANLiMLiOo _ (aS ) 


AMLO% = Miio (40a) 


2 
If Li2O2(s) and Li2Oçs) are simultaneously deposited as the solid discharge product, as described by Eq. (3d), being formed in 
the porous composite cathode; then, the amount of Li2O2(s) and Li2O¢s) deposited for the equimolar formation situation is given, 
respectively, by: 


—AMii a 

AMLO) = (==) MLio, = (45 Mitio: (40b) 

and 
—ANii qg 

AML 06) = (==) Miio = (45 Mitio (40c) 

For the case of different molar formation rates of the solid product species simultaneously: 
1 a 
AMLO) = EF (EF) Mind, (40d) 
AMLin0 4 = —— a Mio (40e) 
ae +r \ aF i 


In the above equations, Mi, Mtipo and Miizo, are the molecular weights of lithium, lithium oxide, and lithium peroxide, respectively 
[ggmol7']. 

For Li2O2s) as the only discharge product, Eq. (3b), the increase in volume of the LizO2,) discharge product in the cathode 
volume is given by: 


AMLO ; max Mii 
AVio = = (4 a (408) 
PLi2016) 2F PLi202%) 


For Li20¢s) and Li202(s) as the discharge products, Eq. (3d), the volume increase of the discharge products in the cathode volume 
for the equimolar formation situation is given by: 


AMLO% P AMLO (E) Miro é Miizo» 
PLi20%) PLi202%) 4F PLizOs)  PLi202%) 


AV(Li206) +Lir02%)) = (40g) 


where PLi,0,) Nd PLiz0.@), are the respective densities of Li2OQis) and Li202¢) [g cm~?]. For the case of different molar formation 
rates of the solid product species simultaneously, the increase in volume of the discharge product species is given by: 


qi" r Miio 1 MLi,0> 
AV (Lin0()+Li202)) = ( T (40h) 
ee ORG 2F l+r PLi206) l+r PLi2O2(5) 


If Odiachs [cm], is the maximum thickness of the solid discharge product layer around every carbon grain present in the porous 
composite cathode reaction layer; then, the volume available for occupation by the solid discharge product around every carbon 


grain is given by: 


4r dave,C- grain max 3 dave,C- grain 3 
Vaischg product on C-grain = 3 2 F Odischg z Zz f dischg product deposited (41a) 
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4r 
2 z ax \2 
= (5 (Orische LO: loaves gain + 1.5dave,C-grain dische F (ôdischg) ] fdischg product deposited (41b) 
Vaischg product on C-grain has the units of [cm?], Jaischg product deposited 18 the factor representing the fraction of a carbon grain surface 
available for the discharge product deposition, and dave,C-grain is the average diameter of the carbon grains present in the cathode. 
The number of carbon grains, each of diameter, dave,C-grain, needed in the cathode reaction layer is given by: 


A Vaischg product(s) (42) 


Nc- grain = 
Vaischg product on C-grain 


Note that A Väischg product(s) = AVL in 0x5) for LigOas) as the only discharge product and that A Vdischg product(s) = AVLix0(6) +Li20245)) 
for Li2O,s) and LizO2()as the discharge products. 
The mass of a carbon grain of diameter, dave,C-grain, 18 given by: 


T 
MC-grain = (=) (Gi ye,C-grain)PC-grain (43) 
where $c-grain is the carbon grain mass density [g m7?]. 


The total mass of carbon grains in the cathode reaction layer is given by: 


total 
™MC_grains = ™C-grain NC-grain (44) 


The total volume of carbon grains in the cathode reaction layer is given by: 


total 
2 C-grains 
total g ( 4 5) 


C-grains = 
oe PC-grain 
The mass of the binder in the cathode reaction layer is given by: 
Mbpinder = Vbinder(Mass of carbon grains in the cathode reaction layer) = Fbinderm aren (46) 


where Ppinder is the ratio of the mass of binder to the mass of the carbon grains in the cathode reaction layer. 
The volume of the binder in the cathode reaction layer is given by: 


Mbinder 
Vbinder = (47) 
Pbinder 
where (binder 1S the binder density [g cm73}. 
The mass of the polymer electrolyte in the cathode reaction layer is given by: 
shit i total 
Mpolymerelectrolyte = Fpolymer electrolyte (mass of carbon grains in the cathode reaction layer) = rpolymer electrolyte!C-grains (48) 


where rpolymer electrolyte 18 the ratio of the mass of polymer electrolyte to that of the carbon grains in the cathode reaction layer. 
The volume of the polymer electrolyte in the cathode reaction layer is given by: 


polymer electrolyte 


Vpolymer electrolyte = (49a) 
Ppolymer electrolyte 
where polymer electrolyte 18 the polymer electrolyte density given by: 
Melectrolyte + 1 
Ppolymer electrolyte = (49b) 


(Melectrolyte/ Pelectrolyte + 1/ polymer) 


where Melectrolyte is the mass of the electrolyte, such as LiPF¢ in an organic solvent mixture, e.g., propylene carbonate and 1,2- 
dimethoxyethane, per unit mass of a polymer, e.g., the polyacrylonitrile (PAN) holding the electrolyte. Ppolymer electrolyte aNd Ppolymer 
are the densities of the polymer electrolyte and the polymer, respectively [gcm7>]. 

It is assumed for this cell construction that some pore volume is available in the cathode for the air or oxygen access to the reaction 
sites even when the cell undergoes maximum discharge. The volume of the pore space available for air in the cathode layer at the 
end of the desired maximum discharge is: 


cathode 
Vair, end discharge = J: air end discharge Viotal (50) 


where fair,end discharge 1$ the volume fraction of the total cathode volume available for the presence of air in the pore space at the end 


of maximum discharge, and where vaio is the total cathode volume [cm?], that must accommodate the desired maximum cell 
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discharge. The total cathode volume is now expressed in terms of the required component volumes as: 


oe = (volume of the solid discharge product(s)) + (volume occupied by carbon grains) 


+ (volume occupied by binder) + (volume occupied by polymer electrolyte) 
+ (pore volume occupied by air at the end of maximum cell discharge) (51a) 


Writing Eq. (51a) in mathematical terms gives: 


yo = A Väischg product(s) + vE“ lains +F Vbinder + Vpolymer electrolyte + J: air,end discharge va (51b) 
where A Väischg product(s) = AVLin025) for Li2Ozs) as the only discharge product when the discharge occurs by Eq. (3b) and 
A Vaischg product(s) = A V(Liz0()+Li2026)) for Li2O,s) and LizO2s) as the discharge products when the discharge occurs, for exam- 
ple, by Eq. (3d) for the situation of equimolar formation of the solid product species. Eq. (51b) leads to the following expression for 
the estimation of the total cathode volume: 
cathode A Vaischg product(s) + Versa + Vbinder + Vpolymer electrolyte 
total = 1 , , (52) 
( _ J air,end discharge) 


One should note that all the terms in the numerator of the right-hand side of the equal sign in Eq. (52) are dependent on the desired 
cell maximum discharge and other parameters of the cathode reaction layer design, including the carbon grain average diameter and 
the maximum thickness of the solid discharge product film formed around each carbon grain. 

The total pore volume fraction occupied by the gas phase in the cathode reaction layer before the start of the cell discharge, p 
is given by: 


total 
fees _ pore A Vaischg product(s) 


pore “ y,cathode ~~ cathode 
V Viotal 


+ fi air,end discharge (53) 


total 


3. Computed results 


(a) Provided data: 
e Concentration of the salt LiPF¢ in an organic solvent such as a mixture of propylene carbonate and 1,2-dimethoxyethane = 1 M. 
e Temperature of cell = 298.15 K. 
e Dry oxygen partial pressure = 0.21 atm. 
(b) Summary of the partial computed data: 
e Reversible cell voltages: 
o Reaction (3b): 3.16 V. 
o Reaction (3d): 3.01 V. 
e Effect of temperature on the reversible cell voltage at constant pressure of 1 atm and Tọ = 298.15 K: 
o Reaction (3b): E= E? — (0.0005)(T — Ty) =3.16 — (0.0005)(T — To) [V]. 
o Reaction (3d): E= E? — (0.0006)(T — To) = 3.01 — (0.0006)(T — To) [V]. 
e Effect of pressure on the reversible cell voltage at constant temperature and gas phase composition (P = absolute pressure, 
atm; T= 298.15 K): 
o Reaction (3b): 


dE RT\ 1 0.0128 
(=) = ( ) = [V atm7!]. 
oP Tyo, 2F) P P 


o Reaction (3d): 


dE 15RT\ 1 0.0096 4 
2 = = [V atm" ]. 
OP) Tyo, 4F P P 


e Equilibrium reaction extent at 25°C and 1 atm: 
o For reaction (3b): Assume ¢/No,,0 = 1, where & is the extent of reaction at the equilibrium state with respect to reaction 
(3b) and No, ,o is the initial g mol of oxygen. Fractional conversion, if oxygen is the limiting reactant (in a batch reaction 
vessel), is: 


XO, 


= No>,0 = Noy,eq =1 No, ,eq =| (=a = =) = Ee 24 (100%) 
No>,0 Nos,0 No3,0 Noo,0 
o For reaction (3d): 
Assume 
e 


o1 
Nono 1.5 
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and 
Noy,0 — Noo,eq 
XO> SS 


1 
No,,0 No>,0 
e Mass depletion and increase rates: 
© igeom = 0.2 MA cMgeom [10]. 
o Total system mass increase rate is then: 
dmt 
sys 
dt 
= 2.49 x 1078 le s7! RA , for reaction (3d). 


o Lithium metal consumption rate is: 


= 1 (100%). 


No», eq = (1 =) z 1.5& 
No,,0 02,0 


= 3.32x 1078 [g s7! cmm]: for reaction (3b). 


—A Mig ~9 —1 am -—2 : 
ie i = 2.07 x 10 [g mols CMgcom]; for reactions (3b) and (3d). 
o Solid discharge product production rates are: 
dNLi O 
-m = 1.04 x 107? [g mol s7! eet for reaction (3b). 
= 5.2 x 107! [z mois! aN , for reaction (3d). 


e Cell heat removal rate for its isothermal operation at 25°C, 1 atm: 
oD = igeom = 0.2 x 107° Acm, 7, [10]. 


geom 


o ES = 2.65 V (a typical value). 
o ỌrxnGb) = —0.00013 and Ọ,xn (Ba) = —0.0001 [J s7! cm; ám]. 

e Theoretical discharge and maximum thermal efficiency: 
o Theoretical discharge per gram of Lis) consumed is 3863.03 mAh g7! Liis). 
o Maximum thermal efficiency, 7thermal,max = 90% for reaction (3b) and =90% for reaction (3d), using the thermodynamic 


data from [11]. 
4. Concluding remarks 
The developed/assembled mathematical formulation has been presented to compute the following: 


e Reversible single electrode and overall cell reaction voltages for the formation of only Li2O2(s) as well as both Li2O¢s) and Li202¢s) 
as the solid discharge product species in the cathode reaction layer. 

Moles of species Lis), LizO,s), and Li2O2(s) as a function of time during the cell discharge period. 

Heat transfer rate between the cell and its surroundings for a desired total current at the actual cell voltage. 

Theoretical discharge capacity and thermal efficiency of the cell. 

The composite carbon cathode volume of the cell. 


Some illustrative cell numerical data have been presented on the reversible cell voltages, the effect of temperature on the reversible 
cell voltage, the effect of pressure on the reversible cell voltage at constant temperature and gas phase composition, the equilibrium 
reaction extent at 25 C and 1 atm, cell mass depletion and increase rates, cell heat removal rate for its isothermal operation at 25 °C, 
1 atm, and the theoretical discharge and maximum thermal efficiency. 

It is envisaged that the developed formulation will be useful in the design/development as well as in the performance analysis of 


a lithium/air or oxygen battery system. 
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